On the attractive plasma-charge system in 2-d 
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Abstract 

We study a positively charged Vlasov-Poisson plasma in which N 
negative point charges are immersed. The attractiveness of the system 
forces us to consider a possibly unbounded plasma density near the 
charges. We prove the existence of a global in time solution, assuming 
a suitable initial distribution of the velocities of the plasma particles. 
Uniqueness remains unsolved. 
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1 Introduction 

An interesting physical situation in Plasma Physics is vi^hen a system of N 
heavy charged particles (say positive ions for instance) evolves in a plasma, 
that is a sea of light particles, of opposite sign (say electrons). The latter 
subsystem is often conveniently described in terms of a mean-field approxi- 
mation by a continuous distribution f{x, v, t), being x, v, t position, velocity 
of a light particle and time respectively. Thus the time evolution of the full 
system is given by the following Vlasov-Poisson equation 

{dt + v-V^ + {E- F)- V,)/(x, V, t) = 0, (1.1) 

where E is the self-consistent electric field generated by the continuous 
charge distribution / (see details below) and F is the electric field generated 
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by the point charges whose positions at time t are denoted by .^^(t), . . . , S,^ (t), 
that is 

AT 

F{x,t) = J2F'{x,t), (1.2) 

i=l 

with 

Here d = 2,3 denotes the dimension of the physical space and we are assum- 
ing, for notational simphcity, that charges and masses of the point charges 
are identical and unitary. 

Equation (|l.ip has to be complemented by the ODE describing the mo- 
tion of the point charges which is, for any i = 1, . . . ,N, 

N 

f,' = -E{e,t)+ F\e,t), (1.4) 

where rf = denotes the velocity of the i-th charge. 

The purpose of this paper is to study the existence of a global solution 
to eq.ns (|l.l|) - (|1.4p . In the completely repulsive case, the issue of global 
existence and uniqueness of the solution has already been approached and 
solved: in [T] and [2] for bounded and unbounded two-dimensional plasma 
distributions and in [11] for a bounded plasma in three dimensions. The 
background on which these papers and the present one are based is the com- 
plete and satisfactory theory for the usual Vlasov-Poisson equation (namely 
equation (jl.ip with F = 0) which has been developed in many articles as 
[SliaElElElliailSllISllISllIT] and others. Of course adjoining point 
charges to a continuous charge distribution implies adding singular forces as 
the field F defined in (|1.2p - ()1.3p . Hence the existing theory is strongly per- 
turbed and has to be deeply modified. Some papers related to this context 
are P El H [IS] . 

The essential tool employed in [1] and [TlJ, which will be used also here, 
is the study of a function /i(x, u, t) defined, in case of a single charge, as 
follows: 

h{x,v,t) = + a G{\x - m\), 

where G is the fundamental solution of the Poisson equation and cr = ±1 
for the repulsive and attractive case respectively. It represents the energy 
of a plasma particle in the reference frame relative to the point charge, 
and of course it is not a time invariant function. Nevertheless it presents 
two essential properties: first, its time derivative does not depend upon the 
singular forces, which is crucial in proving its boundedness, and secondly, 
in the completely repulsive case (cr = +1), it has good sign properties to 
give a control on the velocity of the plasma particle and its distance from 
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the charge. In the model discussed in this paper, that is in the attractive 
case, the second property is not satisfied, being a = —1. Consequently the 
sign of h is not defined and hence getting a bound on h does not imply 
that the single terms appearing in its definition are bounded. Indeed, in 
contrast with the repulsive case, here there can be plasma particles that, 
even starting far apart from the charge, arrive close to it in a finite time 
while gaining arbitrarily large velocities, but still having bounded energy. 

We stress that we succeed here in proving the global existence of the 
time evolution of system p.l|) - (|1.4p . but not its uniqueness, for which we 
think that new considerations are needed. 

The paper consists of six sections. After the introduction, in Section [2] 
we pose the problem and present the main result (Theorem [1]) , stated for a 
system consisting of a positive plasma density and a single, negative point 
charge. We start with the single-charge case in order to provide a proof 
which is clear and contains all essential tools that are needed also for the 
A'^-charges case. Sections [3] to [5] are then devoted to the proof of Theorem 
[11 In Section [3] we introduce a family of regularized differential systems, 
for which we establish many preliminary estimates holding uniformly with 
respect to the regularization. In Section H] we show the main technical result 
in this paper, that is the boundedness of the above mentioned function h 
(see Theorem [2]). This allows, in Section [5l to prove the existence of a 
global solution of the system, obtained as limit of the regularized dynamics. 
Finally in Section [6] we state and prove global existence of a solution to 
system (jl.ip - (jl.4p for N charges (Theorem [3]) . 

2 The result for a single charge 

In this section we consider a plasma in two dimensions with only one charge. 
We set i(,{t),r]{t)) for position and velocity of the charge at time t, being 
(^, rj) their initial data. Moreover f = f{t) e L°°{L^ r]L°°) denotes the den- 
sity of the plasma, and we assume that /(O) = /o is a bounded probability 
density. 

Equation (jl.ip describes a conservation law for the density along the 
time evolution of the characteristics, which is, at least formally, given by 
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the following differential system: 



x{x, V, t) = v{x, V, t) 

v{x, V, t) = {E - F) (x(x, V, t),t) 

{x{x, V, 0),v{x, V, 0)) = (x, G \ {^} X 

f — y 

E{x,t) = _ p{y,t)dy 



Jm2 \x-y 




< 



(2.1) 



p{x,t) = / f{x,v,t)dv 
Jr2 




f (x(x, V, t),v{x, V, t),t) = fo{x, v), 



together with the evolution of the charge, moving according to: 



'm = vit) 
< fi{t) = -E{m,t) 

^(e(0),7?(0)) = (C,7?)GM2xM2. 



(2.2) 



Note that, if /o is smooth, then any solution {^{t),r]{t); f{t)) to system (|2.ip - 
([2:2]) satisfies the Vlasov-Poisson equation pTT]) - (fLll) . Clearly, the ODE in 
()2.ip are not well-defined if some plasma particles collide with the charge in 
finite time; However Theorem [1] below ensures that, despite the attractive 
interaction between plasma and charge, such collapses can be essentially 
avoided under suitable assumptions on the support of /q. 

We introduce the function: 



which represents the energy of a plasma particle in the reference frame of 
the moving charge. Moreover we set Sq for the support of /o, which can be 
possibly an unbounded set, and we define the quantity 



with C a sufficiently large constant for further purposes. We will prove 
that if ^(0) is finite, then T-L{t) remains finite on bounded time intervals 
and consequently the velocities of the plasma particles are logarithmically 
diverging as they approach the charge. Nevertheless such slight divergence 
will not prevent us to prove global existence of a solution to (|2.ip - (j2.2p . 



h{x,v,t) = - 



V -r]{t)\^ + In \x -^{t)\ 



n{t) = sup sup \h{x{t),v{t),t)\ +C 



se[o,t] ix,v)£So 
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In the sequel we will often use the notation 

ln_ r = — Inr x(r G (0, 1]), 

with x(^) the characteristic function of the set A. 

We will set C for a positive constant and Q, i = 1,2,..., for some 
constants to be quoted in the course of the paper. All of them will possibly 
depend on ||/o||l°°, ||/o||li and on an arbitrarily fixed time T. Finally, for 
sake of brevity we will sometimes use the shortened notation {x{t),v{t)) 
instead of {x{x,v,t),v{x,v,t)). 

Our main result is the following: 

Theorem 1. Let (^,77) G ^ m2 ^ L°°(m2 x M^) be a probability 

density supported on the set 

So = {{x,v) G X . \h{x,v,0)\ < Co} (2.3) 

for some positive Cq. 

Let T > 0. Then there exist 

f eL"^ ([0,r];L°° nL^(M2 X R^)) , E G ([0, T]; (M^)) , 
(C(-),r/(-))GCi([0,T])2, 

and for duo-i-i- ix,v) G Sq there exists {x(-),v{-)) G C^([0,T])^ such that 
(^x{t),v{t);^{t),r]{t); f{t)) satisfy system (|2.ip - (j2.2p on [0,T]. Ln particular, 
for djiQ-a.a. {x,v) G So and Vt G [0,T] we have \x{t) — > and 

f{x{t),v{t),t) = fo{x,v). (2.4) 

Moreover, 

n{T) < C (2.5) 

and finally 

\p{x,t)\<C{l + ln.\x-Cit)\), (x,t) G X [0,r]. (2.6) 

Remark 1. In the proof of Theorem we shall use the fact that (|2.6p . 
which is a consequence of (|2.5p . ensures the uniqueness of the solutions 
{x{t),v{t);^{t),rj{t)) to the ODE in ^^-^2^) once f{t) (hence E{t)) is 
given (see Corollary and Lemma However notice that we are not 
claiming the uniqueness of the triple [x{t),v{t);^[t),r]{t); f{t)^ . 

Remark 2. We stress a significant difference with respect to the repulsive 
case treated in ^T] and J 11']- In those papers the assumption for 71(0) to be fi- 
nite was equivalent to assuming a finite distance between charge and plasma 
at time t = 0, while in the present case it is not so. This is an intrinsic dif- 
ficulty in this setup since, even assuming an initial positive distance between 
plasma and charge, we could not exclude that some plasma particle arrive 
at any prefixed distance from the charge. 
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Remark 3. Hypothesis ()2.3p on the support of Jq ensures that its spatial 
support is bounded. Indeed we have |x — < e'^"~2'^~''' < e*^". By ()2.5p 
this property will be preserved in time (see (|4.23p ). 

Remark 4. We do not claim that the bound (j2.6p is optimal. 



3 The approximating system 

In this section, we introduce a regular version of the original system (|2.ip - 
(|2.2p by mollifying the singular field F created by the charge. More precisely, 
for a small parameter < e < 1, we consider the smooth increasing function 
lue : [0, +oo) — ;> M such that 



ln=. r > 2 In e if r < e. 



lUf r = In r if r > e. 



(3.1) 



We consider next the unique solution ((x£(t), f£(t); ^(.(t), ^(t); /^(t)) to 
the following e-problem on [0, T] : 



Ve = {Ee - Fe){Xe,t) 
{Xe{0),Vem = {X,V) 

= -E,{ie,t) 
(e.(0),r/,(0)) = (e,r/) 



(3.2) 



where: 



Ee{x,t) 
Fe{x,t) 

Pe{y,t) 



Pe{y,t)-r^ ^dy 

\x — 

Fe{x,t) = V^ln,\x-Ut)\ 



fe{y,v,t)dv 



and satisfying: 



fe{Xe{t),Veit),t) = foix,v), 



(3.3) 



(3.4) 



where /q is a smooth, compactly supported approximation of /q. Here we 
are in presence of the Vlasov-Poisson problem with an additional smooth 
gradient external field, for which the classical theory for global existence and 
uniqueness of the solution applies with minor modifications. 

Thanks to ()3.4p we have ||/£(i)||/^oo = ||/o||l°o. Moreover, since p.2p is 
hamiltonian the flow (x, v) i— )• (^^(x, v, t), Ve{x, v, t)) preserves the Lebesgue's 
measure on x M? for all t > 0. Consequently all the ||/e(t)||LP norms 
{p > 1) are conserved. 

We introduce next the regularized relative energy per plasma particle: 



he{x,V,t) = hv - rie{t)\^ + lUs \x - ^s{t)\ 



(3.5) 
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and we set: 



Tisit) = sup sup \he{Xe{t),Vs{t),t)\ + Ci. (3.6) 

se[o,t] ix,v)eSo 



Notice at this point that definition (j3.6p does not allow us to consider initial 
data /o satisfying assumption (|2.3p , since there are configurations of particles 
for which h{x,v,0) is bounded while /ie(x,f,0) is not. This is due to the 
smoothed potential In^ which cannot compensate large velocities of particles 
that are very close to the charge. To overcome this difficulty we introduce 
another positive parameter /3 > e and, instead of considering initial data /g, 
we consider /q supported in the set 

= {{x,v) G m2 X 1^2 . \h{x,v,0)\ < Co,|x-e| > /?} C So. (3.7) 
We observe that Sq is a bounded set (see Remark [3]). Setting now 

^f(t)= sup sup \he{Xe{t),Ve{t),t)\+Ci (3.8) 

we have by definition (j3.7p : 

n^{0) = n{0) < Co (3.9) 

and /q is a compactly supported function satisfying the assumptions of 
Theorem [TJ 

We will prove that the solution fs{t) to system (|3.2p with initial condi- 
tion /q enjoys estimates independent of e and j3 making it possible to pass 
to the limit as e — )• and /? — )• 0. To simplify the notation from now on 
we will sometimes omit the index /3, but we keep in mind that the solution 
depends on both parameters. We emphasize that all constants appearing in 
what follows do not depend on e and (3. 

Let us introduce the total energy of the system, which is a conserved 
quantity: 



^"^l^ - y\Peix,t)pe{y,t)dxdy + j \ns\x - ie{t)\Pe{x,t) dx. 



(3.10) 



Our first observation is that <?e(0) is positive and bounded uniformly in 
£ and /3, as it is stated in the following 

Proposition 1. If fo is supported on the set given by ()3.7p . then every 
single term appearing in the definition o/i5e(0) is bounded. As a consequence 
we have: 

< £e{0) < Cs 

provided that C2 has been chosen large enough. 
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Proof. We start by recalling the following elementary fact: \/p > there 
exists a positive constant C{p) such that 

j {In-rfdr = C{p). (3.11) 

Next, for any (j;, v) G Sq we infer from definitions ()3.5p and ()3.ip that 
<2(|t;-r?|2 + |??|2) 
< 4(/ie(x, V, 0) + ln_ \x - ^1) + 2|??|2. 

Therefore we deduce from definitions (|3.8p - (|3.9p that, provided the constant 
Ci is sufficiently large, 

\v\'^ <C{n{<d)+\n.\x-i\) V(x,v)e5^. (3.12) 

This implies that for x G we have 



Pq{x) = I fo{x,v)dv 

<\\M\o.J x{\v\ < Cy^nO) +ln_\x - ^\) dv, 

hence 

Po{x) <C{'H{0) + ln.\x-^\), VxGM^. (3.13) 

Estimates (|3.12p and (|3.13p enable us to prove that any single term 
appearing in the definition ()3.10p of £{0) is bounded. Indeed, by Remark [3] 
on the compactness of the support of po and by (j3.13p and (j3.1ip we have 

J I lue \x - (,\\poix) dx < C. (3.14) 

Moreover, for the same reason: 



<c(^l + Jln_\x- ^1 (ln_ Ix - CI + 1) dx^ 



hence 

\v\'^ fo{x , v) dx dv < C. 



Finally: 

ln|x - y||po(2;)/0o(y) dx dy 
<C I ln|x - 2/||(ln_ \x - C\ + l)(ln- \y - C\ + l)dxdy. 

J J supppo 
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Again by ()3.1ip the above integral can be easily bounded by means of 
Cauchy-Schwarz inequality, so that <f^e(0) is bounded uniformly in e and 
positive, provided C2 is sufficiently large. □ 

The preceding result does not give us any e-uniform bound on the single 
terms composing since it could be bounded uniformly in e by com- 

pensation. The next two results provide such informations. Their proof is 
extensively based on the conservation of the Lebesgue's measure and on the 
invariance of the plasma density along the motion of the characteristics. 

We set 

Ke{t) = \ I \v\^Mx,V,t)dxdv+^-^^^. 

Proposition 2. 

sup Ke{t) < Ca, (3.15) 
te[o,T] 

sup / \x\p^{x^t) dx < C^, (3.16) 
te[o,ri J 



sup ||p,(t)||i2 <C76. (3.17) 

t6[0,T] 



Proof. For any M > 0, we have 



Pe{x,t)= I fe{x,V,t)dv+ / fe{x,V,t)dv 

l\v\<M J\v\>M 



<TrM^\\f,{t)\\L^ + ^ J \v\^f,{x,v,t) 
By optimizing in M we find: 



dv. 



1/2 

Pe{x,t)<C{ I \v\'^fe{x,V,t)dv 



whence, by definition of Kf,{t): 

\\pe{t)\\L^ < C/Kjt). (3.18) 

On the other side, from the energy conservation and Proposition [1] it 
follows that 

Keit)<C3 + j \n^\x-^sit)\psix,t)dx 

+ 7;[[ ^ri-\x - y\pe{x,t)ps{y,t)dxdy. 
^ JJ\^-y\>i 

Now, Cauchy-Schwarz inequality and (|3.18|) yield 

sup y ln_ \x - z\ps{x,t)dx < C\\ps{t)\\L2 < C^K^(t). (3.20) 



(3.19) 
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Moreover < Inr < r for any r > 1. Hence (|3.19p and ()3.20p imply: 

(3.21) 



K,it)<Cs + C^/K;it) + - 1 1 {\x\ + \y\)pe{x,t)pe{y,t)dxdy 



< C3 + C^KS) + j \x\pe{x,t) 



dx. 



Next, by (|3.4p and the fact that the flow preserves the Lebesgue's mea- 
sure on X we have 

\x\p£{x,t) dx = I |xe(x, w, t)|/o(x, f ) dx 

X (3.22) 



< / I + / \v£{x,v, s)\ds \ fQ{x,v) dx dv. 







Recalling Remark[3l po has compact support, so that applying again Cauchy- 
Schwarz inequality we get 



(3.23) 



j \x\pi;{x,t) dx < C + j j \v\fe{x,v, s) dx dv ds 

<c[i + j'^^/KA^d^ . 

Going back to ([321]), (13:23]) implies: 

K,{t) <C + C^KS) + C f VKe{t)ds 

Jo 

and ()3.15p follows then from Gronwall's Lemma. Finally recalling ()3.23p 
and (IXTBD we conclude that UiSAM and (IXTTP fohow from (f3l3|) . □ 

Remark 5. We stress the fact that ()3.15p yields a bound on the velocity of 
the charge and consequently on its motion, which remains confined over the 
interval [0, T]. 

Proposition [2] implies the following bounds on the potential terms in the 
energy: 

Corollary 1. 

sup \lni;\x - ^£{t)\\pe{x,t)dx <C 
te[o,T] J 

sup // \ ln\x -y\\pe{x,t)pi;{y,t)dxdy <C. 
tG[0,T]J J 
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Proof. Arguing as in ()3.19p - ()3.23p . by Proposition[2]we have thanks to (|3.20p 
and ([3:23]) 



/ K \x - ^e{t)\Pe{x,t)dx 

J\^-Ut)\<l 

+ - 11 lii\x-y\pe{x,t)ps{y,t)dxdy <C. 

^ J J\x-n\>\ 



'\^-y\l 

On the other side from definition (j3.10p of the energy it follows that: 



\\ie\x - ie{t)\pe{x,t)dx ^ \ 11 In^ \x - y\ps{x,t)pe{y,t) dx dy 



< C3 - / \ne\x - S,e{t)\pe{x,t)dx 

J\x~eJt)\<i 



+ -// ln\x - y\peix,t)peiy,t)dxdy. 

'\x-y\>l 



2 



(3.25) 

Hence the conclusion follows from ()3.24p and ()3.25p . □ 



4 The function 

The main result of this section is the following 
Theorem 2. 

ne{T) < c. 

The proof of Theorem [2] requires some preliminary results, which are 
stated hereafter. 

Lemma 1. For all {x,v) G and for any t € [0,T] it holds: 



\Ve{t)\ < C^Ueit) + \n.\x,{t) - Ut)l (4.1) 
Pe{x,t)<C{ne{t)+\^-\x-Ut)\), (4.2) 



and 



m < c^JnM. (4.3) 



Proof. We infer from Remark [5] that ()4.ip and ()4.2p are the equivalent to 
(|3T2D and (f3l^ at time t. 

Now we prove (|4.3p . Recalling Remark [3] in Section [51 let B{0,R) be 
the ball of radius i? > 1, so large that it contains the support of po and so 
that moreover ^^(t) G B{0,R),yt £ [0,T]. Let us fix t G [0,r] and define 
t* = max{s G [0,t] : Xe{s) G B{0,2R)}. Then either t* = t ov t* < t and 
in the latter case we have |xe(s) — ^sis)\ > ii > 1 for any s G {t*,t]. Hence 
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in view of the definition (|3.5p of he it follows that for s £ {t* , t] we have 
< y^2'H£{s), which implies 

\Xe{s)\ <2R + j^^ \Ve{T)\dT < C^JUeis) 

provided the constant Ci in (j3.6p has been chosen sufficiently large. □ 
Proposition 3. For t € [0, T] we have 



ll^eWIlL- <C^/^n7^tj■ 
Proof. We decompose Ei;{x,t) as 

Es{x,t) =Ii{x,t) +Z2{x,t), 

where 



f — y 

Ii{x,t)= Pe{y,t)- 12 

J\x-v\<s F-yr 



l2{x,t) = pe{y,t)-^ — ^dy 

J\x~y\>5 k-yr 

and < (5 < 1 is to be determined hereafter. Let us first estimate the term 
Ii{x,t). We have by (lOl : 

|2i(x,t)|<C / . . dy 

J\x-y\<5 \x-y\ 

<Cneit)6 + c[ '""-ly-^f^ dy. 

J\x-y\<5 \x-y\ 

By considering the two cases: 

\y - Ce{t)\ < \x - y\ and \y - > \x - y\ 

we arrive at 

f ^^-\y-^^(')\ dy< [ ^-^dz<C6\ln6\, (4.4) 

J\x-y\<S \x-y\ J\z\<S \z\ 

SO that 

\Iiix,t)\< Cne{t)6 + C5\ln6\ <Cne{t)6 + C. (4.5) 

On the other side, applying Cauchy-Schwarz inequality and using (j3.17p we 
get: 

P6{y,t) , , f Pe{y,t) 



\U^,t)\<r 'f^dy+[ 

J5<\x-y\<l F ~ y\ J\x 



I i^y 

I5<\x-y\<l \X — y\ J\x-y\>l F — 2/1 



<\\PemL^[ ^ +1<Cv4M. 

Vi5<|2|<l Pl / 



12 



With the choice 

we obtain the thesis, provided the constant Ci in the definition ()3.8p is large 
enough. □ 

Proposition 4. For t E [0, T] we have 

sup / \w\ ^'^y^'"^f> dydw < Cy^ne{t)ln'He{t). 

xeR2 J \x -y\ 

Proof. Let < (5 < 1 to be chosen hereafter. We perform the integral as 
follows: 



11 fe{y, )J_ _ 2^^ri.^ f\ _|_ X2(x, t) + Is{x, t), 

F - y\ 



with 



i{x,t) = [ 

J\x 



fe{y,w,t ) 
\x-y\<s ■ ■ \x-y\ 



Xi{x,t) = / \w\ — — ^ — - — dydw 



-r ( .\ I' II fe{y,w,t) 
l2[x,t) = / \w\ — i — dydw 

Js<\x~y\<l F ~ y\ 

T I f II feiy^'w^'t) , , 

i.3(x,ij= / \w\ — j j — dydw. 



\x-y\>l 



\x - y\ 



Recalling estimates ()4.ip and ()4.2p and proceeding analogously to Proposi- 
tion El we get: 

J\x-y\<5 \x-y\ 



<c5neitf/' + c [ dz 

J\z\<S \z\ 



Next, using Cauchy-Schwarz inequality, estimate p.lSp on the kinetic energy 
and ()4.2p we obtain: 

1/2 / , , ,x \ 



fe{y,w,t) 

S<\x~y\<l 



\x — 



M"^ fe{y,w,t)dydw] \^ I 7.! 7 dydw 

\ 1/2 

Pe{y,t ) 
s<\x-y\<i \x — y\ 



<c{ I T^^dyj 

1/2 



\J5<\x-y\<i k-yr 



\ 1/2 

dy] 
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Now, arguing as in in ()4.4p we have 

f , ln_ \y-^e{t)\ ^ 

j<5< x-v <i F ~ yi 



'<5<|x-y|<l 

hence 



X2(x,t) < C V^^MM + C"! lnJ|. 
Finally, using again ()3.15p we arrive at 

T-i{x,t) < C. 

Hence we conclude that 

I \w\ /^(^'^'^) ^y^^ < C6'He{tf''^ + CJ'He{t)\\^5\+C\\n5\ + C 
J \x-y\ 

and the thesis is achieved by choosing 5 = l-Leii)^^ ■ D 

Now we state a " quasi-Lipschitz" property for the field E^, which is a 
modification of a standard inequality (see, e.g., [lOj). 

Proposition 5. We have for t G [0, T] and x, y G 

\E,{x,t) - E,{y,t)\ < Cifilx - y\){ne{t) + ln_ \x - y\) (4.6) 

where 

ip{r) = r(ln_ r + 1). 
Proof. Let d = \x — y\. If(i>l/3we apply Proposition [3] to obtain 



\E,{x,t)-E,{y,t)\ <2\\E,{t)\\L^ < C^JlnUeit) < Cd^JlnUeit). (4.7) 
Otherwise we set z = {x + y)/2 and we make the following decomposition: 
\Ee{x,t)- Ee{y,t)\ < Ii{x,y;t) + l2{x,y;t) + h{x,y;t), 

where 

Ii{x,y;t)= [ (^^73^ + ^^) 



h{x,y;t) 



\z-z\<2d 



2d<\z~z\<^ 



1 1 



\x — z\ \y — z\ 



P£{z,t) dz 
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By (|4.2p . proceeding as in ()4.5p we get 

h{x,y;t)<C [ {ne{t)+ In. \z- urn dz 

J\z-x\<3d F ~ ^\ 

+ C I (neit) + In. \z-Ut)\) dz 

J\z-y\<3d \y - Z\ 

r 111 
< Cd'He{t) + C / ^^^-Y^dz, 

thus 

h{x,y]t)<Cd{ne{t)+\n^d). (4.8) 
For the second integral we write, always by ()4.2p : 

Pe{z,t) 



h{x,y;t) <Cd [ 
J2, 



I 1 2 

2d<\z-z\<^ 1^ ~ ^1 



dz 



i2d<|2-2|<i k ~ ^1 

< Cd-Hsit) In^ d + Cd [ ln_ \z-Ut)\ 

J2d<\z^z\<^ \Z-Z\'^ 

We further split the integral above into two parts, obtaining on the one 
side 

Li -„i^^^^^l#^^<ln-d / —^dz<C{ln^df, 

J2d<\z-z\<- \z-z\^ J2d<\z--z\<h\z- 



\z-Ut)\>d 

and on the other side 



L, ^"-1^ dz<^ f In^ \z -Ut)\dz<C In. d 

J2d<\z-z\<\ \z--z\2 - M'^ Ju-iAt)\<d ' ^'^ ~ 



|2-Ce(t)|<d 

Hence, gathering both estimates we get: 

I \n.\z-Ut)\ ^ ^ ^) ^ 

whence 

l2{x,y;t) <C{ns{t)d\n.d + dlnld). (4.9) 

Finally for the term I3 we observe that if |z — z| > 1/d then minjlx — 
z\,\y- z\}> l/{2d). Hence: 

^3(2;, y, t) <Cd j pe{z, t) dz = Cd. (4.10) 

Estimates (03]), (gj]) and (fiJOjl imply the thesis. 

□ 
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Before presenting the proof of Theorem [2] we state a useful identity: 

Lemma 2. 

t 

Vs{s) ■ Es{Xeis),s) ds 

= <^>s{xe{t),t) - 0) + /" / My, w,s)w ^j""^ ~ ^ dy dw ds, 

Jo J \Xe{s)-yr 

where $e is the potential due to the plasma, that is 







<^six, = y In |x - y\ pe{y, t) dy. 



Proof. Since Es{x,t) = Vx^eix,t) we have: 

d d 

Ve{t) ■ Ee{Xe{t),t) = - — t), 

which imphes: 

/■* /■* d 

I Ve{s)- Ee{Xe{s),s)ds = ^e{Xe{t),t)-^{x,{))- I —^e{Xe{s),s)ds. 

(4.11) 

We next evaluate the partial derivative in ()4.1ip . For any z G M^, we have 
d d f 

—^e{z,s) = — fe{y,w,s)ln\z -y\dydw 
d f 

TT h{y,w)\n\z -ye{s)\dydw 



ds 



Z - Vejs ) 
\Z - ye{s) 



fo{y,w) Weis) ■ — — -^j-^dydw. 



Hence the thesis follows. □ 

Now we are in position to present the 
Proof of Theorem 2. 

Let {x,v) S Sq and consider the characteristic (^^(t), ^^(t)) starting at 
time t = from {x,v) and its relative energy hi;{x£{t) ,Ve{t) ,t) for t G [0, T]. 
We compute 

jh,{Xe{t),vS).t) = (Veit) - Ve{t)) ■ t) + ^^(6(0, *)) (4-12) 

and we see that the singular part disappears in the derivative. We next look 
for an estimate for the time derivative of /i^ in terms of quantities regarding 
the plasma which we have all already estimated. 
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Let 6{t) be a continuous function, to be chosen later, such that < 
6{t) < 1. Assume that at time t € [0, T] we have \xe{t) -Ce{t)\ > S{t). Then, 
recalhng that supfgp'p] < (see ()3.15p ). we infer from Proposition 

[3] and from (jH]) that 



^he{Xe{t),Ve{t),t) 



<2\vS) - r^e{t)\\\ES)\\L^ 



(4.13) 



< CyJUe{t) + ln_ 5{t)yJ\n'He{t)- 

If on the contrary \xe{t) — CE(i)| < (5(t) then we write eqn. ()4.12p as: 

jh,{Xe{t),vS).t) = {Ve{t) - 7]e{t)y{E,{Ut),t) - E,{Xe{t),t)) 

+2{Ve{t)-r]e{t))-E,{Xe{t),t). 

We start by estimating the first term in ()4.14p . for which we can use the 
quasi-Lipschitz property stated in Proposition [5] that is: 

\ve{t)-r], (t) 1 1 E, (e. (t) ,t)-E,{Xeit),t)\ 



< C^Ue (t) + ln_ I X, (t) -Ut)\\Ee{Ut),t)-E, {Xe {t),t)\ 

< C {Ueit) + ln_ \xe{t) - Ut)\?'^ - Um 

< c-He{tf^^{\x,{t) - um + c (in_ \xe{t) - ut)\f' vii^eit) - um- 

Now, by definition of the function ip we have y?(r) < ip{6{t)) for < r < S{t), 
and moreover 

(^(r)(ln_rf < C Vp > 0. 

Therefore we obtain 

\ve{t) -rj,{t)\\E,{Ut),t) - E,{x,{t),t)\ < C {ne{tf'\{m) + l) • (4.15) 
Hence from (|il3]) . (jilil) and (fiJ5]) it fohows that: 

\h{x,{t),Ve{t),t)\ 

<niO) + C [ yJUeis) + 5{s)yJ\nUe{s) ds 
Jo 

+ C j (ne{sf'\{5{s)) + l) ds + 2 I {v,{s)-r],{s))-E,{xe{s),s) ds 

(4.16) 

The last integral on the right-hand side can be estimated by p.lSp . Propo- 
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sition [3] and the preceding Lemma [21 obtaining: 



{Ve{s) - r]e{s)) ■ Ee{Xe{s),s) ds 



< 



Veis) ■ Es{Xeis),s) ds 



+ C [ ^JlnU^is) ds 
Jo 



<C [' ^Inneis) ds+\^e{Xe{t),t)-^{x, 0) 

Jo 



+ 



fs{y, w,s) w ■ j^^^ ^ dy dw ds 



\Xe{s) - 



This by Proposition S] imphes: 



{Ve{s) - rjs{s)) ■ Ee{Xe{s), s) ds 



< \^e{Xe{t),t)\ + mx,0)\+C f C^-He{s)\n-He{s) ds. 

Jo 



(4.17) 



It remains to estimate the potential terms in (|4.17|) . Arguing as in (|3.14|) . 
we infer from Remark [3l ()3.13p . ()3.1ip and Cauchy-Schwarz inequahty that 
|$(x,0)| < C. Moreover by ([3:20]) and Proposition [2) 



(4.18) 



Pe{y, t) ln„ \xe{t) -y\ dy < C. 
On the other side estimates ([3.16P and ([4.3p give: 

/ pe{y,t)\n\xe{t)-y\dy< Pe{y,t){\xe{t)\ + \y\) dy 

<C + \x,{t)\ I pe{y,t) dy<Cy^7Ut), 
which, together with estimate ([4.18p . prove that 



\<^eiXe{t),t)\ < C (7^4(0 + ij < Cy/UAt)- 

The use of ([TO]) in (lilTD yields 



(4.19) 



{ve{s) - r]e{s)) ■ Es{xs{s), s) ds 
<C^Jne{t)+C I Cy^'He{s)ln'He{s) ds. 

Jo 



(4.20) 



Now, inserting ([4.20p in (|4.16p we have for any {x,v) S S^: 



\hix,it),veit),t)\< n{o) + c^/7Ut) 

+ C [ yJUeis) + ln_ 5{s)yJ\n'He{s) ds 
Jo 

+ C f (ne{s f^ip{S{s)) + y''He{s)lnns{s)) ds. 
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Finally, taking the supremum over {x, v) G Sq we conclude that 



+ C [ v/ln^e(s)ln_ 6{s) ds 
Jo 

+ c (ne{sf'^^{5{s)) + ^/ne{s)\nn,{s)) ds. 

By choosing 
we are led to: 

nsit)<nio) + c^neit) 

+ C {luUeis) +ne{sflHnne{s) + ^H,{s)hlHe{s)) ds ^^^l) 

<n{0) + Cy/HjJ) + C [ -Heis) ds. 

Jo 

If the constant Ci in the definition ()3.6p of 'He(t) is large enough, we finally 
obtain: 



ne{t) < Cn{0) + C [ -Heis) ds (4.22) 

^0 



with the constants not depending on e and /3 and the conclusion follows 

from Gronwall's Lemma. □ 

Setting Sf for the support of the density /^(t), the following corollary is 
a direct consequence of Theorem [2j 

Corollary 2. There exist positive constants independent of e and (3 for 
which it holds: 

Sf = {{x,v) E X m2 : l^l < C, ne{t) < C} (4.23) 

Pe{x,t)<C{l + ln^\x-Ut)\)- (4.24) 

Moreover: 

sup ||^e(t)||Loo < C (4.25) 

te[o,T] 

and for any {x, y) G Sq : 

sup \Ee{x,t)-Eeiy,t)\<C-fi\x-y\) (4.26) 
te[o,r] 

where 

7(r) = r(2 + ln_ r)^ if < r < 1; j{r) = 4 otherwise. (4.27) 
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Proof. All the bounds are direct consequences of Tlieorem[2l Estimate ()4.23p 
follows from (|iT3]) in Lemma [H (|4.24p and (|4.25p follow from (fOj) and 
Proposition [3] respectively. Finally ()4.26p follows from Proposition [5] in case 
1^; — y| ^ 1 and from (j4.25p otherwise. □ 



Remark 6. Notice that the function 7 is positive, increasing, continuous 
and concave on M^. 

5 Proof of Theorem [1] 

In this section we prove the convergence of the regularized system introduced 
in Section HI To do this we keep first /3 fixed and state the e-convergence 
results in the following Propositions [6] and [71 Then, using the fact that all 
our estimates are uniform in e and f3 we will be able to remove also the 
cutoff /3. 

We consider the solution to system 

(|3.2p - (|3.4p for some sequence e„ — t- as n — J- cxd. The following result holds. 

Proposition 6. There exists a subsequence of {en} (which we still denote 
by {en}) and there exists E G C(M^ x [0,T]) such that Efr^ converges to E 
uniformly on the compact sets o/M^ x [0,T]. 

Proof. By ()4.25p the sequence is uniformly bounded on the compact sets 
of X [0,T]. Moreover, in view of ()4.26p . it is uniformly equicontinuous in 
X. The result will be a consequence of the Ascoli-Arzela's theorem, once we 
have proven that it is also uniformly equicontinuous with respect to time. 
To this aim we choose t, s G [0,T] such that |t — s| < 1 and we introduce a 
positive, bounded increasing C°°-function g defined as: 

g{r) = 1 if r > a, g(r) =0 if r < — , (5.1) 

where a G (0, 1) will be suitably chosen hereafter, with the further property 
that for some C > 1/2 

< g{r) < -. (5.2) 
a 

Next, we write 

E,^{x,t)-E,Jx,s) 

[1 - g{\x - y|)](pe„(y,i) - Pe,Xv^ s))- ^ dy 

\x - 

9i\x-y\){pe„{y,t) - PsAVis))-. T^dy 

\x - 

= Ii{x;t,s) + l2ix;t,s). 
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By definition of g and ()4.24p we have 
\h[x; t,s)\ < sup ' 



dy < C sup 

&[0,T] J \x-y\<a \X - y\ t&\Q,T] J \x-y\<a 



1 + ln, \y-ieAr) 
\x - y\ 



Hence by (j4.4p we obtain: 



sup sup t, s)\ <C [a + a\ \na\) < Ca\ lna|. 

Q<t,s<Tx&? 



Now we estimate the term l2{x;t,s), writing 



h{x]t,s) = j fo{y,w) 
Since 

d x-ye„{T) 

dr |x-y£„(r)|2 

we have 



d f n r \\\ ^~ VsAt) 

^('"l"-^'-'"'l' |.-.,„(r)P 



dT 



(5.3) 



dy dw. 
(5.4) 



"-^"^ +^M^-^^4;^(x-..„(r)), 



\x-yeAT) 
d x-ye^T) 



dr |x -?/e„(T)|2 



\x-yeSr)\^ 

k£n(^)l 



Hence from ()5.4p and properties (|5.ip - ()5.2p for g( it fohows that 



\h{x;t,s)\ < J foiy,w) 

C f 
+ — / fo{y,w) 



dr 



dy dw 



< 



C 

c 

a 



foiy,w 



' \weAr)\ 
* \^e,Ar)\ 



\w\ 



X-ye„{T 
fen{y,W,T 



■dT 
-dT 



dy dw 
dy dw 



dy dw dT. 



\x - y\ 

Thanks to Proposition [3] and Theorem [2] we conclude that 

\l2{x;t,s)\ < C 



(5.5) 



and the conclusion follows from (|5.3p and ()5.5p by choosing a = \t — s\^^'^ . □ 

Proposition [6] is an important step to prove the uniform in time conver- 
gence of (^£„ , ) and of {x^^ {x, v) , Vs„ {x, v)) for almost-every initial config- 
uration {x,v). However we need also to control the size of the "bad initial 
configurations", those leading to possible collapses with the charge at some 
time. This is done in Lemma [3] below, the proof of which is postponed in 
the Appendix at the end of this section. 



dy. 
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with C independent of n and /3. 

Finally, the following statement is the last result we need in order to 
prove Theorem [U 

Proposition 7. The sequences x^^ andvs„ are Cauchy sequences in {dfj.o;C{[0,T])) . 
The sequences and rj;,^^ are Cauchy sequences in C([0,T]). 

Proof. Let us set 

Xn,m{t) = " Xe^{t)\ + l^eji) " ier.{t)\. 

We have for all t G [0, T] : 



Fix now a positive parameter 6 such that 2max(e„,em) < 5 < 1. We 
decompose Sq as 



with SniS) the set introduced in Lemma [3l Clearly, for fixed S the fields 
-Fe„ are bounded and Lipschitz uniformly with respect to n on T{6y; This 
make it possible to handle the last integral in the right-hand side in ()5.6p . 
On the other side. Lemma [3] provides a control on the size of the bad set 
T(5). Therefore in order to prove Proposition [7| we will first let n, m — )■ oo 
for fixed 5, then let 5 — )• 0. 
To that aim we introduce 




(5.6) 



s^ = mur{6r 



where 



r{6) = Sn{6)USm{6) 
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and we notice that by Lemma[3]/io(r((^)'^) — t- 1 as 5 — t- 0. 

We start by estimating X^^{t). We estimate the first two terms in the 
right-hand side of (|5.6|) by observing that 

\EeAXeAr),r)-Ee^{Xe^{T),T)\ 

< I ixe„ (r) , r) - E,^ (r ) , r) | + | (r) , r) - E,^ (r ) , r) | 

(5.7) 

and by making an analogous decomposition for the second term. By the 
quasi-Lipschitz property ()4.26p and Remark [6l we have 



\Ee^ixeAr),r)-Ee^{xe^{T),T)\ + \E,^iCeAr),r) - Ee^iC,M,r) 



(5.8) 



Next, if (x,f) belongs to T{5Y, we have on [0,t] 

Xsn (t) - S,e„ (r) (r) - (t) 



|F,„(xe„(r),T)-F,„(x,„(T),r)| 
therefore 



ke„(r)-C,„(T)|2 \xUr)-^eArr 



\FeAXe„{T),T)-F,^ix,jT),T)\ < 



|^n,m(''") 

6^ 



(5.9) 



Going back to (j5.6p . in view of (j5.7|) . (j5.8p and (|5.9p we obtain for all 

{x,v) G r{sy 



<C sup \\E,St)-E,^{t)\\L^+C f r 
te\o,T] Jo Jo 



^n,m(''") 
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(ir ds 



C 

< C sup \\Ee,M - E,^{t)\\L^ + -2 / / 7(^n,m(r))drds 

t6[0,T] Jo Jo 

SO that, being 7 an increasing function, it follows 

xl^{t)<c sup ||ii;,„(t)-ii;,„(t)||icx> 



+ 



te[o,T] 

c /■* 




7o Jo 



/r(<5)= \ne[o,T] 



7 sup Xn,rn{u)\ d^lQ{x,v) 



dr ds. 



By concavity of 7 we may apply the Jensen inequality to the above integral, 
obtaining: 



Xl^{t)<C sup \\E,St)-E,^{t)\\L^ + ^ [ f\{Xl^{T))dTds. 
te[o,T] Jo Jo 

(5.10) 
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Next, we have in view of Proposition [6l 

lim sup \\E,„{t) - E,^{t)\\L^ = 0. (5.11) 
"'™^°^te[o,T] 

This, by definition of 7, imphes that the second order integral inequahty 
()5.10p can be handled by choosing n, m large in function of 5. More precisely, 
defining 

uj{n,m)= sup \\Es„{t) - Es^{t)\\L-- 
te[o,T] 

we obtain 

<^(t)<C[..(n,m)]-P(-^*/^). (5.12) 

The proof of ()5.12p is elementary and is given in Lemma H] in the Appendix 
at the end of this section. We want to stress that it is possible to prove it 
only because (|5.1Up is a second order integral inequality. 

Let us now estimate the remaining term X^^{t). In view of the bound 
(|4.23p on the support of the spatial density of the plasma and of Lemma [3] 
we get: 

Xl„,{t) < Cfio (r(5)) < C6\ \nSf/\ (5.13) 
Gathering ()5.12p and ()5.13p we obtain 

Kmit)+<mit) < C[a;(n,m)]-P(-^*/^) +C<5|lnd|3/2. (5.14) 

Hence, defining 

we conclude that 

limsup;f„,m(i) < C6\lnSf/'^. 

n,m—>oo 

This, by arbitrariness of 6, implies the first part of the thesis, that is 

is a Cauchy sequence in (d/zo; C{[0,T])) and is a Cauchy sequence in 

c([o,r]). 

For the velocities we introduce analogous definitions, that is: 

"^«.™(*) = / '^^P ^n,m{s) dfloix,v) 

^l,mii) = / sup Vn,mis) dHoix,v) 

Jr{5) s€[o,t] 



and 
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Proceeding in analogy with the previous computation, we infer from (|5.12p 
that 



Vl^it) < Cu{n,m) + ^ ^ 7«™(r)) dr 

< Cuj{n,m) + ^7([w(n,m)]^"P(-^*/^)). 



(5.15) 



On the other hand, Cauchy-Schwarz inequahty combined with the bound 
on the charge velocity yields: 

(t) 

/ \ 1/2 

<C|r(<5)|l/2 f sup {\VeAs)\' + \Ve^is)\'') foix,v)dxdv] +C|r(5)| 

= C\T{d)\'/^( [ \v\^ sup if,„{x,v,s) + fe^{x,v,s)) dxdv] +C\T{d)\. 

\J se[0,t] J 

Therefore, thanks to the bound (|3.15p on the kinetic energy and to Lemma 
El again we obtain 

Vl^{t)<C6'/'\ln6\'/\ (5.16) 
Again, in view of (|5.1ip . ()5.15p and ()5.16p we are led to 

limsupV„,m(t) < C(5i/2|ln5|3/4 



n,m—^oo 



and the conclusion follows as before. Hence the proof of Proposition [7] is 
complete. □ 



The results achieved up to now allow us to complete the 
Proof of TheoremUi 

Thanks to Proposition [3 there exists {S,i')iV{')) ^^^^ ^^i dfiQ -a. a. (x, v'j , 
there exists {x{x,v, ■),v{x,v, •)) such that (■?£„(•); ^en(")) converges to {C{-),t]{-)) 
and, for dfiQ-a.a. (x, v), (x£^(x, v, •), t'e„ (x, v, •)) converges to {x{x, v, ■),v{x, v, •)) 
uniformly on [0,r]. It follows that the map {x,v) i— t- {x(x,v,t),v{x,v,t)) 
preserves the Lebesgue's measure on Sq for all t G [0,T]. 

Next, for t S [0, T] we define the measure diJ,{t) = {x{-, •, t), v{-, 'i i))^ dfio 
on M2 X R"^, i.e. 

j ip{x,v)dn{t) = j ip{x{t),v{t)) fo{x,v)dxdv V99 G ^(M^ xR^). (5.17) 

Clearly, by (|3.4p and (|5.17p the sequence fe„{t) converges to d^{t) in 
the weak sense of measures for all t S [0, T]. On the other hand, fe„{t) is 
transported by a measure-preserving flow and it satisfies the uniform bounds 
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(|3.15p and ()4.23p : therefore, it is also uniformly equi-integrable and Dunford- 
Pettis theorem ensures that /e„(t) is weakly relatively compact in L^. As a 
result, we have dfi{t) = f{t)dxdv for some f{t) S L^. Since, on the other 
hand, /e„ is uniformly bounded in ([0,T],L^ n L°°) we conclude that 
/ G ([0,T],Li nL°°). 

We may now define p = J f dv G L°°[L^). By the same arguments, we 
check that Pe^it) converges weakly in to p{t) for t G [0, T]; in particular 
the bound ()2.6p in Theorem [1] holds for p. Setting then E = p* x/\x\'^ and 
using the fact that p satisfies (|2.6p . we obtain, mimicking (for example) the 
arguments of the proof of Proposition [3l that E'e^ converges to E uniformly 
on X [0,r]. 

Furthermore, Lemma [3] shows that the set {{x,v) G Sq : inf^gjo.T] l^i^) — 
£,{t)\ = 0} has zero d^o-^sasure. This ensures uniform convergence of the 
singular field for dfiQ-a.a. initial data. We can then pass to the 

limit in equations (|3.2p to find that {x{t),v(t);^^{t),r]{t);f(t)) satisfy eqns. 
()2.ip and ()2.2p over [0,T]. In particular, (|2.4p follows from ()5.17p and from 
the fact that {x,v) i— )■ {x{x,v,t),v{x,v,t)) preserves Lebesgue's measure. 

As a matter of fact, we proved Theorem [1] only for /q being supported 
in Sq. However, for any /o supported in Sq, let us introduce the sequence 
of initial conditions 

f^{x,v)=Afpfo{x,v)x{\x-^\> ^) 

where JV/b is a normalization factor since we are working with probability 
distribution. Obviously A/}? — )• 1 as /3 — )• 0. Let f^{x, v, t) be a corresponding 
sequence of solutions we have already constructed. Since all our estimates 
are uniform in /3, we can remove the /3 cut-off exactly the same way as we 
just did with e, achieving thereby the proof of Theorem [TJ 

Actually we could also have proven our result in a different way, working 
with a single sequence, by choosing /3 = /3(e) suitably vanishing with e. 

Appendix 

Proof of Lemma [3] 

For < 5 < 1/4 we set Sn{6) = S^^) U Sl{d), where 

Sli6) = {ix,v) e Sni6) : X e B{^,25)} 
Sl{6) = {{x,v) G S^{6) : X € B{^,26y}. 

The estimate on the measure of S^{6) is trivially given by the hypothesis 
(j3.7p on the initial data and its consequence 

\Sl{S)\< [ x{S^) dx dv < C [ ln_ \x - ^\ dx<C 5^ \ln5\. 

(5.18) 
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On the other hand, let (x, v) G S^{6). By continuity there exists to = io{x, v) 
such that \xs^{to) — Ce„(*o)| = ^- We set G [0,T] for the connected 

component containing Iq such that 5/2 < \xs„(t) — ^e„(t)| < 25 for t G 
(t~,t+). By virtue of (|4.ip and Theorem [2] for any {x,v) G S^{5) we then 
have 

\veAt)\^ <C\ln5\, VtG(t-,t+), 

therefore 

\t+ -t-\ > r ■= AT. 

Now we partition the interval [0, T] into N{5) + 1 intervals [tj, tj+i] of length 
smaller than AT/2. Then {fjt'^) has to contain at least one of the ti, so 
that 

N{5) 

Sl{5) C U '5n,i(5), 
i=0 

where 

SnM = {{x,v) G : 5/2 < k.Jti) - ^eAti)\ < 2<5}. 

Next, since the flow {x,v) i— )■ (xe„(x,f,tj),fe„(x,f,tj)) preserves the Lebesgue's 
measure on x we have 

\SnA^)\ = \{{x,v)GSn{t.i):5/2<\x-^e,Ati)\<25}\ 

where is the support of fs„{'^i)- Thus, in analogy with estimate (jS.lSp . 

by Theorem [2] we get: 

\Sn,i{5)\ < C5^\ln5\. 
Since by construction N{5)AT < CT, we finally obtain 

N(5) 

\Sli5)\ < ^\Sn,ii5)\ <CN{5)5'^\ln5\ <C5\\n5\^/^, (5.19) 

1=0 

and the conclusion follows from ()5.18p and ()5.19p . 

□ 

Lemma 4. Let u be a positive, continuous function on [0, T] such that 

u{t) <h + a I I 7(ti(r)) drds 
Jo Jo 

where 7 is the function defined in (|4.27p . a > and < b < 1 such that 

b < exp(l -e^^^). 
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Then 

u{t) < v{t) 

where the function v{t) is the solution to the differential equation 

ij = a j{v) 

v{0) = b v{0) = 0. 

Moreover: 

v{t) < C76exp(-2Tv/5)_ (5 20) 

Proof. The proof of the first part of the statement is standard so that we 
are left with the proof of ()5.20p . By the properties of the function 7 there 
exists T* <T maximal such that 6 < v{t) < 1 on (0,T*). By multiplying 
both terms in the differential equation by i) we get for t G (0, T*) 

2dt^ ' dt^ ' 

where (/> is a primitive of 7. Notice that, by the definition of 7, (j) is an 
increasing function such that (l){v) < 2f^(l — Inu)^ Mv G [0, 1]. Therefore 



V < V2a(0(?^)^^^(6)) < ^/2a(p{v) < 2^f^ v{\ - Inu). 

Since the function ^(1 — Inu) is positive and increasing in (0, 1] we can apply 
the Gronwall lemma to this differential inequality, getting: 

v[t) < exp(l - e-2rv^)6^-P(-22^VS)^ t E (0,r*). 

By the choice of b and the definition of T* , we obtain T* = T, so that the 
previous inequality holds on [0, T] and the Lemma is proved. □ 



6 The case of A^-charges 

This section is devoted to the system already presented in the introduc- 
tion, consisting of negative point charges and a positive plasma. Setting 
r7*(t)) for position and velocity of the i-th charge at time t with initial 
condition (^*,?/*), the equations (|l.ip - (|1.4p in terms of characteristics are: 

'x(t) = v{t) 

v{t) = E{x{t),t) -F{x{t),t) 
x(0) = X, v{0) = V 

< ^^(*) = ^^(*) (6.1) 
r,^it) = -E{e{t),t)+Y^F^e{t),t) 

r(o) = r, v(o) = v, z = i,...,iv 

J{x{t),v{t),t) = fo{x,v) 
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where E, F and have been defined in ()1.2p and (|1.3p . 

In order to present our global existence result we need some notations 
to describe the support of the initial density /o in system ()6.ip . Let do be 
the minimal distance between two charges at time t = 0. We set 



4 



N 



i=l 



where B{C, do/4) = {x : \x - f | < 4/4}. 

We introduce the energy of a plasma particle relative to the i-th charge: 

h\x,v,t) = ^\v-r^\t)f +\n\x-e{t)\. 



We further define 

ho{x,v) -- 

We assume that the support of /q is the set given by 



\h\x,v,0)\ if {x,v) G Aj, 
max 0)1 if {x,v) £ A'^. 



So = So (co, y) = {(^'^) ■ ho{x,v) < Co} (6.2) 
for some positive Cq. 

Remark 7. In order to have a finite Cq in definition ()6.2p we need to 
decompose the set of initial data by means of the sets Aj and A'^. This is 
due to the fact that whenever a particle is close to the i-th charge, then its 
velocity has to be large because h^{x,v,0) is assumed to be bounded. On the 
other hand, for the same particle, h^{x,v,0) would diverge for any j ^ i, 
since it is far from the j-th charge and its potential part cannot compensate 
such a large velocity. 

We are now in position to state the main result of this section: 

Theorem 3. Let . . . be distinct points of M."^, r]^ , . . . ,rj^ G and 
fo G L°°(M^ X M^) be a probability density supported on the set So{Co,do/A) 
for some positive Cq. Let T > 0. Then there exists 

f eL"^ ([0,r];L°° nLi(M2 X R^)) , ^ G L~ {[0,T]; L°° (R^)) , 

ieit-),v'i-))eC\[0,T]f, i = l,...,N, 

ix{-),v{-)) G C\[0,T]f fordfio-a.a. {x,v) G So{Co,do/A) 

such that (2;(t), f (t); ^*(t), T/*(t); /(t)) satisfies system ()6.ip on [0,T]. In par- 
ticular, for dfio-a.a. {x,v) G 5*0 and Vt G [0,T] we have \x{t) — ^*(t)| > 
for all i and 

f{x{t),v{t),t) = fo{x,v). (6.3) 
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Finally, 



< C(l +maxln_ -f (t)|), (x, t) G x [0, T]. (6.4) 

i 

To prove Theorem [3] we establish some a priori estimates assuming that 
a smooth solution to system (|6.ip does exist. This will allow us to complete 
the proof by means of a regularization procedure as the one we have used 
in Section [5] in the case of a single charge. 

The total energy of this system reads now 

^(i) = 2 / \v\''f{x,v,t)dxdv + Y.-^l^- ^ ln\eit)-e{t)\ 

i=i i^j=i 

N 

(6.5) 



1=1 



This is a conserved quantity along the solutions to system ()6.ip and, as 
before (see Proposition [T|), the assumption (j6.2p on the support of /o ensures 
that £{0) is bounded and positive for a choice of the constant C large enough. 
However we cannot a priori exclude that, by compensation, a couple of 
charges moves closer while another couple, or the plasma, go to infinity. 
Actually our strategy for the proof of Theorem [3] consists first of all in 
establishing a lower bound for the mutual distances between the charges. 
This will imply that a plasma particle can approach at most one charge at 
the same time and consequently the A'^-charges case will be reduced to a 
sequence of one-charge problems, which we already solved. In order to show 
that the N point charges remain separated during the motion over the time 
interval [0, T] we introduce the kinetic energy of the system 



N 



K{t) = ^ f \v\'^fix, V, t) dxdv + '^ 

i=i 

and we prove the following result: 
Proposition 8. 

sup K{t) < C5. 
te[o,r] 



2 
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Proof. The bound can be proven by paraphrasing Proposition [2j Indeed, 



K{t)=£{0)+ ^ ln\e{t)-e{t)\ + - [ ln\x -y\p{x,t)p{y,t)dxdy 

N 



N 

<£iO) + ( V yin\e{t)-eit)\ + l: [[ ln\x -y\p{x,t)piy,t)dxdy 
^i^l ^JJ\x~y\>l 

N 

+ / In- k-f (t)|p(x,i)cix 

(6.6) 



where X] j stands for the sum over ah i ^ j such that \^l{t) — > 1. 

The last two terms in (j6.6p can be estimated as we did in (j3.19p in Propo- 
sition [2l so that by ()3.23|) we get: 

1 ^ 

- II In \x - y\p{x,t)p{y,t)dxdy + y~] / \n.\x - C{t)\p{x,t) dx 

^JJ\x-y\>l f^J 

<C\\p{t)\\L2 + l- [[i\x\ + \y\)p{x,t)p{y,t)dxdy 



< CyjK{t) + j \x\pe{x,t)dx 

< C^/Kif) + C I ^/K{s)ds. 

Jo 



(6.7) 



Besides, we have 

N N 

yin\e{t)-e{t)\<Y{\e{t)\ + \em 

i^j=l i,j=l 

N 



<C+iZ I ds{W{s)\ + W{s)\) (6-^ 

<C+ [ yjK{s)ds. 
Jo 



Finahy, combining ()6.6p . ()6.7p and ()6.8p we obtain: 

r-t 



K{t) < £{0) + C(l + 7K(t)) + ^ y^K{s)ds. (6.9) 
We conclude by means of Gronwall's Lemma. □ 
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As a consequence of Proposition [8] we have the result we were looking 

for. 

Corollary 3. There exists d = d{T) > such that: 

min inf ICHt) - ^Ht)\ > d. (6.10) 

i^j t&[0,T] 

Proof. By definition ()6.5p of the energy we have: 

N 

- ln\e{t)-e{t)\ = -C, + £{0)-K{t) 

N 



1 

+ 2 



/ / - y\pix,t)p{y,t)dxdy - / ln\x - Cit)\p{x,t) dx 



< £"(0) + ^ // ln\x -y\p{x,t)p{y,t)dxdy 

'\x-y\>l 



2 

N 



+ ^^-\x - C{t)\p{x,t)dx. 

1=1 •' 

(6.11) 

But \x — y\ > 1 implies < ln|x — y| < |x| + |y|, so that by (|3.23|) and 
Proposition [8] the first integral in (|6.1ip is bounded by a constant. For the 
second one we recall the bounds (|3.2U|) . (|3.18|) and again Proposition [5] to 
conclude: 

'"ln_ \x-e{t)\p{x,t)dx < C. 



Hence we have proved that 

N 

- \^\e{t)-e{t)\<c. (6.12) 

Now notice that the bound on the kinetic energy implies 

sup \e{t)\<c, 2 = i,...,iv, 

tG[0,T] 



SO that, putting as before X] j to indicate the sum over all i ^ j for which 
If (t) - i^{t)\ > 1, by (f6l^ we have: 

N N N 

Y in^\e{t)-e{t)\ = - Y inie(t)-e^(t)i + ( ^ )*inie(t)-c^'(t)i 

i^j=l ii^i=^ i¥'j=^ 

N 

< ( Y yin\e{t)-e{t)\ + c 

N 

<2Y\e{t)\ + c<c 
1=1 
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and this implies the thesis. □ 

Corollary [3] is the key to pass from one to N charges and to establish the 
Proof of Theorem O 

We denote by St the support of the density f{t). Next, we set : 
Kit) = B [e{tl -J X K^t) = n Af(t) 

where d is the constant in (|6.1U|) . and 

%i{t) = sup sup \K'{x,v,s)\ + C 

se[0,t] (x,i;)GAi(s)n5s 

T-L^it) = sup max sup f , s)| + C, 

se[0,i] * (a;,j;)eA=(s)nSs 

with C sufficiently large for further purposes. Finally, let 
nit) = max i max-?^i(t), rL\t) 



Remark 8. Since do/4 > d/S, there exists Ci > such that 5o(Co,do/4) C 
5o(Ci,(i/8). Hence ()6.2p implies thatT-L(Q) is finite. 

Clearly, we have 

AT 

\v\'^ <2'H{t) + C^\n_\x-e{t)\, VtG[0,T], V(x,v)GSi; 

In particular, one can readily transpose to the present case the a priori 
bounds given, in the single-charge case, for the electric field in terms of the 
largest energy (Propositions O [Hand [SJ. Hence, as before (Theorem [2D, our 
goal is to prove that 1-L{T) is bounded in term of the initial datum. The 
conclusion of Theorem [3] will then follow by simply mimicking the arguments 
of Section [SJ we omit the details here. The strategy in obtaining an a priori 
bound on T-L{T) is the following. We consider a small time interval [0, A*] 
given by 

A* = / (6.13) 

t* < T being a suitably small fixed time (see ([6])). We will prove that %{t*) 
is bounded, by making an iteration procedure in time over the intervals 
[{k - 1)A*,A;A*],A; = l,...,n = {[t*] + 1)/A*. Since the time t* is only 
function of constants, the subsequent prolongation of the estimate from t* 
to T will be straightforward. 
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Now we claim: 



Claim: During the time interval [0,A*] any characteristic of the plasma 
starting from 5o(Co,do/4) may approach at most one charge. More pre- 
cisely, if there exists an index i for which 

x{so) G B ( f (so), ^ ) for some sq £ [0, A*], (6.14) 



then 

xit)eB(eit),^] VtG[0,A*]. (6.15) 



Indeed, (l6J5]) clearly holds if 3;(s) G B (f (s),(i/8) on [0, A*]. Otherwise 
we have |x(si) — C('5i)| = d/8 for some si. We denote by (s_,s+) the 
maximal connected component of si in [0, A*] on which 

^ < |x(s) -f (s)| < y, sG(s_,s+). 
Let s £ (s_,s+). If \x{s) -f (s)| > d/8 then x{s) G A%s), so that 

\v{s)\ < y/wJF) + c < yjm{t*). 

Otherwise, we have x(s) G Aj(s) but \x{s) — ^*(s)| > d/lQ hence we obtain 
again 

\v{s)\ < ^2\h'{x{s),v{s), s)\ + C^l In \x{s) - e{s)\ < ^/wJF). 
By choice of A* (see (j6.13p ) we thus obtain 

\\x{s) - e{s)\ - ix(si) - e{si)\\ < VmF)A* = ^ 

lb 

We conclude that (s_,s+) = [0, A*] and (fOSD follows. 

Now let us consider t G (0, A*], an index i and (x, v) G Aj(t) n St- Then 
it has to be 

{x,v) = {x{x,v,t),v{x,v,t)) 

for some {x,v) G Aj(0) U A^(0) since we just proved that a characteristic 
may visit at most one ball i?(^*,d/8) on [0, A*]. 
For s G [0,t] we compute 

^^h\x{s),v{s),s) = ivis)-v\s)) ■ {Eixis),s)+Eieis),s)) 

N 

-{v{s)-rj\s))- J2 {FHe{s),s)+F^{x{s),s)). 
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Notice that, as in the case of one charge, the singular contribution due to 
the nearest i-th charge to x disappears. For j ^ i, it follows from (j6.15p 
that 

\F^{e{ts),s)\ + \F^{x{s),s)\<^ 

so that the external fields are smooth gradient fields on [0, A*] and the 
results stated in the preceding sections hold with minor modifications. Con- 
sequently we arrive at the analogue of estimate ()4.2ip : 

\h\x{t),v{t),t)\ <'H{0) + C^n{t) + C [ nis) ds. 

Jo 

Taking the supremum we obtain 

sup \h\x,v,t)\ <n{0) + C^/H(t) + CA*'H{A*). (6.16) 

(x,€)eAi(t)n5t 

Consider now (x, G A'^(t) n St- There exists {x,v) such that 

{x,v) = {x{x,v,t),v{x,v,t)) 

and {x,v) G Aj(0) U A'^(O) for some i. Proceeding as above we get again: 

\h'{x,v,t)\ <n{Q) + Cy/Wt) + CA*'H{A*). (6.17) 

Now we observe that since (x, -u) G A^(t), for any j ^ i we, have 

\h^{x,v,t)\ < \K\x,v,t)\+C\v\ + |ln|x-f (t)| - In |x - ^■''(t)! | 
< \h\x,v,t)\+C^/n(f) 

so that by (j6.17p we find 

max sup \h\x,v,t)\ < niO) + Cy^n{t) + CA*n{A*). (6.18) 
1=1,. ..,N (x,v)eA'={t)nSt 

Hence (j6.16p . ()6.18p and the monotonicity of T-L^t) in t imply for any t G 
[0,A*]: 

'H{A*)< n{0) + Cy^niA*) + CA*'H{A*) <'H{0) + C^/n{t*) + CA*n{t*) 

provided the constant in the definition of 1-L{t) is large enough. Recalling 
the definition ()6.13p of A* we arrive at 

^(A*) < niO) + CA*'H{t*). 

By iterating the procedure n = {[t*] + 1)/A* times we get: 

nit*) < n{nA*) < ^(0) + CnA*n{t*) = n{o) + ct*n{t*) 
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which imphes 

■H{t*) < n{o) ^ 



1 - Ct* 



Hence by fixing 

t* = 1/2C 

we get a bound on Hit*) : 

n{t*) < 2'H{0). (6.19) 

To go from t* to 2t* we iterate the procedure, by dividing the interval [t*, 2t*] 
into n intervals of width A* = d/{16^/3n{2t*)) with n = {[t*] + 1)/A* and 
we find 

^(2f ) < n{f) + ct*n{2t*) 

from which by (j6.19p it follows 

1«(2f ) < 2^(0) = 4^(0)- 

By the choice of t* the iteration stops at Nt* = T so that for = 2C([T] + 1) 
we finally obtain 

niT) < 2^n{o) < c^n{o). (6.20) 

The conclusion of Theorem [3] now follows from ()6.20p and a suitable adap- 
tation of the arguments presented in Section O as was explained at the 
beginning of this section. □ 
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